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Improved actions in SU(2) and SU(3) lattice gauge theories are investigated with an emphasis 
on asymptotic scaling. A new scheme for tadpole improvement is proposed. The standard but 
heuristic tadpole improvement emerges from a mean field approximation from the new approach. 
Scaling is investigated by means of the large distance static quark potential. Both, the generic and 
the new tadpole scheme yield significant improvements on asymptotic scaling when compared with 
loop improved actions. A study of the rotational symmetry breaking terms, however, reveals that 
only the new improvement scheme efficiently eliminates the leading irrelevant term from the action. 
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INTRODUCTION: 



At the beginning of lattice gauge theories, Wilson pointed 
out that it is important to maintain exact gauge invari- 
ance for finite lattice spacings a thus enforcing gauge in- 
variance in the critical limit of the lattice model. The 
minimal choice for an action which satisfies this crite- 
rion is nowadays known as the Wilson action [lj. Re- 
lying on the concept of universality, any lattice action 
which incorporates the correct symmetries in the con- 
tinuum limit should work in principle. However, some 
lattice actions do give better approximations for contin- 
uum results for coarser lattices. This issue is central 
for computer simulations and has led to a continuous 
development of the so-called improved and perfect ac- 
tions H i II i, i, 0, S ©, M M M 0, G3, EE M- 

The basic idea is to add terms to the action which are 
irrelevant in the continuum limit, but which give better 
approximations at finite lattice spacing 0, H| . Different 
proposals have been made for such actions on the basis 
of perturbation theory 14, la |q| or using renormalisation 

MM 



group techniques [11|, [12|, [13 



It is widely ac- 
cepted that the so-called tadpole improvement is impor- 
tant for good properties of these actions @, [l(| • To our 
knowledge, however, a systematic study of different im- 
plementations of tadpole improvement has not yet been 
carried out. 



In the context of computer simulations, an extrapola- 
tion of data to the limit of vanishing lattice spacing is 
necessary. Such extrapolations can be made very trust- 
worthy if a relation to an ab initio continuum calculation 
can be established. If we use the string tension a as 
the fundamental energy scale of Yang-Mills theory, the 
perturbative treatment of continuum SU(N C ) Yang-Mills 
theory predicts the dependence of the lattice spacing on 



(3 = 2N c /g 2 (with g the bare gauge coupling) to be: 
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There, the 1-loop and 2-loop coefficients 
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are universal. The dimensionless parameter c a depends 
on the observable and must be determined by non- 
perturbative methods such as lattice simulations: any 
physical mass scale (call it m) in units of the reference 
scale is independent of the lattice spacing for sufficiently 
large (3 and is obtained from 
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Let N denote the number of lattice points in one di- 
rection of the lattice. In actual lattice simulations, (3 
cannot be chosen too large if we wish to work with rea- 
sonable lattice sizes, Na(f3). Using the Wilson lattice 
action, it turns out that these values of (3 are still too 
small to observe perturbative scaling: for moderate (3 
values, large corrections to the scaling function (JXJ) are 
observed. Note, however, that the ratio m 2 a 2 /era 2 is al- 
most independent of the lattice spacing at these (3 values 
which let us reliably calculate low energy observables. 
This property, called 'scaling' in the literature, must not 
be confused with "asymptotic scaling" , i.e. perturbative 
scaling, which is the focal point of the present paper. 

In this paper, a new tadpole improved action is proposed. 
The construction of this action offers a new understand- 
ing of the otherwise heuristic "derivation" of the standard 
tadpole action commonly used in simulations nowadays. 
We will find that tadpole improvement is highly impor- 
tant for the approach to asymptotic scaling for reason- 
ably sized lattices. Finally, a thorough study of rota- 
tional symmetry breaking effects obtained from the static 
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quark potential will reveal that the standard tadpole ac- 
tion does not cancel the leading order irrelevant terms 
of the action. The numerical data suggest that a com- 
plete cancellation might be achieved by means of our new 
action. 



II. ACTION AND IMPROVEMENT 



A. Standard tadpole improvement 

Because of the non-linear relation between the link field 
Un{x) and the continuum gauge potential, lattice per- 
turbation theory suffers from large tadpole contributions 
which, however, must cancel for an extrapolation to the 
continuum limit. That these tadpole contributions are 
indeed large can be easily anticipated from the expecta- 
tion value of the link in Landau gauge, Uj}(x). A naive 
expansion with respect to the lattice spacing, i.e., 



trUj}(x) 
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implies that this expectation value should be of order 1. 
Actual simulations show, however, that the latter expec- 
tation value is at most of order 0.8 although era 2 is as 
small as 0.05. In order to improve the approach to the 
continuum limit, one considers the ratio between the link 
and its expectation value. One assumes that the devia- 
tion from unity now provides a better definition of the 
(continuum) gauge field A^x): 

^W/^KW) = 1 + iA^x)a ... . (5) 

Since Landau gauge fixing is problematic because of the 
Gribov problem, an ad hoc description for tadpole im- 
provement has become standard: defining 



Ufj.(x) = U^(x)/u 



(6) 



where uq is the fourth root of the plaquette expectation 
value, each link field in a lattice operator should be re- 
placed by Ufj,(x). Note that this procedure is heuristic, 
and that many other choices for tadpole improvement ex- 
ist: one could also choose for uq the 8th root of the expec- 
tation value of the 2x2 planar Wilson loop. Nevertheless, 
the prescription outlined above has become standard. 



B. Motivation of the new action 

Let us consider a quadratic (planar) Wilson loop of side 
length s with an orientation specified by /i, v located at 
site x. A naive expansion of this operator yields (see 



e.g. (2) of 0): 



— Re tr 



(x) = 



Oi{x) 
O e (x) 



1 Re tr 

21° 6S6 + -- 
(D^ v ) 2 (x) - 



(7) 



(8) 



{D v F^ v f{x) 



The subscripts ^i, v at O^q have been suppressed. Al- 
ready terms of order 0(s 6 ) break of rotational symmetry. 
Choosing the minimal length s = a, W coincides with the 
(minimal) plaquette, which is the only term in the Wilson 
action. 

As outlined in the previous subsection, numerically the 
expansion (|7|) converges badly for reasonable lattice sizes. 
Without any simulation, this fact can be also understood 
from continuum perturbation theory: although mani- 
festly gauge invariant, in lattice regularisation the high 
energy modes are cutoff at a momentum scale h-uv ~ 
iv /a. It is well known that in cutoff regularisations ex- 
pectation values such as O4 and Oq in ((SJ) diverge with 
the cutoff: 
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The origin of these divergences are quantum fluctuation 
of the order of the cutoff scale. Obviously, these fluctu- 
ations invalidate the expansion ([7]) (choose s = a for the 
moment). However, they do not spoil the calculation of 
physical observables well below the cutoff scale (as will be 
detailed below). Note, however, that if the desired goal is 
to match with asymptotic scaling provided by continuum 
Yang-Mills theory, an expansion such as (J7J) should cover 
high energy modes too. 

One choice for such an action is obtained by replacing all 
operators O n of the action by 



Or, 



(Or, 



(9) 



Only the deviation of the gauge invariant operator from 
its (potentially) divergent expectation value contributes 
to the action. In practice, this construction can be re- 
alized by considering the ratio between the Wilson loop 
and its expectation value. Using ([7]), one can show to all 
orders that in this case the operators O n only appear in 
the combination ©. 

Without resorting to the naive expansion anymore, we 
now assume that the above ratio has a sensible expansion 
with respect to s: 

RetrW M „(x) / (Re tr W^x)) = 1 - i Re tr 



Oa s 4 
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Note that the term O4 gives rise to the continuum action 
proportional to F 2 . The subject of improvement is to 
remove terms of higher order in s from the action. 
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A popular choice (Symanzik improvement) is to use a 
rectangular 1x2 loop. This scheme invokes yet another 
expansion, i.e., of the rectangular loop with respect to 
s similar to the one in (|10p . and relies on a matching 
of the expansion coefficients to eliminate the irrelevant 
terms. Here, we are going to use a 2 x 2 quadratic loop 
which is a scale transform of the plaquette. The moti- 
vation for this choice is that we need to invoke only one 
type of expansion evaluated at two values for the scale 
parameter. The hope is that, because of the scale rela- 
tion between the lxl and 2x2 loops, the cancellation of 
irrelevant terms is more complete at finite values of the 
lattice spacing where the naive Taylor expansion (such 
as (HU])) appears unjustified. 



Assuming that the expansion p0[) is valid at least for 
s < 2a, we use the expansion for s = a (plaquette) and 
s = 2a (2 x 2 Wilson loop) to get rid of the irrelevant 
terms. Defining 



Re tr W^\x) / (Re tr W^\x)) , (11) 



Re tr W$\x) / (Re tr W$\x)) 
we choose for the action 



S = p 



[«i^)+«s^(4 (12) 



Using the expansion (JTOJ) for s = a and s = 2a, we are 
led to 



16 K 2 
64 k 2 



(13) 
(14) 



This first line ensures compatibility with continuum 
Yang-Mills theory whereas the choice of the second line 
eliminates the order a 6 terms. The solution of the latter 
set of equations is given by 



«i 
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(15) 



The present improvement scheme eliminates from the ac- 
tion contributions from tadpole loops. The main purpose 
for this elimination is that these loops are absent in the ab 
initio continuum formulation of Yang-Mills theory. The 
impact of these loops is therefore to spoil proper scaling 
which is familiar from continuum perturbation theory. I 
point out that, once the tadpole contribution were elim- 
inated, further improvements might be achieved by in- 
voking the standard perturbative improvement scheme. 
We leave such an investigation to future work. Here, 
we will justify by numerical calculations that the new 
action (without further perturbative improvements) al- 
ready gives rise to much better scaling properties. 



the case of standard tadpole improvement, the rule Q 
implies that 

P^x) = ^-RetvW^ix) /u 4 

= Re tr W^(x) / (Re tr W^\x)) , 
= ±-RetrW^ 2 (x)/4. 

While for our new action the numerical burden is a self- 
consistent calculation of the expectations values 

(Re tr W^ l {x)) , (Re tr W$\x)) , 

standard tadpole improvement appears as an approxima- 
tion to this numerical problem: There, only 

(RetrW^ xl (*)> 

is calculated sclf-consistently, and the expectation value 
of the 2x2 Wilson loop is obtained with the help of the 
mean field approximation (in Landau gauge): 



— RetvW^(x)) « u» = K) 



N, 



RetvW^ix) 



Having identified the standard approach as an approx- 
imation to the present scheme, the crucial question is 
whether the properties of our action fully justify the 
higher level of numerical sophistication. The remaining 
two sections will answer this question. 



III. NUMERICAL SIMULATION SETUP 

A. Thermalisation 

The dynamical degrees of freedom are the SU (N c ) ma- 
trices Ufi{x) which are associated with the links of a A^ 4 
cubic lattice. The partition function is given by 



S[U] (wu, W22) 



VU^ exp^K^u,^)} (16) 



..^ L3iuii(/3) F 
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where the action S depends on the expectation values: 



C. Comparison with the standard tadpole 
improved action 

Let us assume that we are dealing with an action which 
features the plaquette and the 2x2 Wilson loop. In 



wnW) = (Retr^ xl (.x)) , 
w 2 2(P) = (RetrW^ 2 (x)) . 



(18) 



Expanding the expectation values of the latter equations 
in terms of their functional integrals, we arrive at a set 
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FIG. 1: Convergence of the action parameters tun and W22 as a function of the number of Monte-Carlo sweeps: SU(2) (left) 
and SU(3) (right). 



of two non-linear equations which determine the two un- 
known parameters wn(0) and 1022(0)' 

w u (P) = \\ W.Retr^ 1 ^) 

exp{S[U}(wu,w 22 )} , (19) 
w 22 (0) = \\ VU^RctrW^ix) 

eMS[U](w n ,w 22 )} . (20) 

Before we can start to accumulate statistically indepen- 
dent lattice configurations {U^} for each 0, we must solve 
the latter set of equations for u>u(0) and w 22 ((3), and we 
must generate a "statistically important" configuration 
by means of thermalisation. 

While the reader is invited to develop their own method- 
ology for this task, we here briefly outline our approach 
which serves the purpose. It appears to be quite natu- 
ral to solve the set of equations (|19I20[) and to generate 
the thermalized configuration within one process. We 
here used a simple iterative procedure: denoting w^(0), 
i = 1, 2, by the approximate solutions to Wu(f3) of the nth 
iteration, better approximations are generated by 

< +1) (/3) = 1 JvU^RetrW^Hx) 

e*p{S[U](w$> t wW)}, (21) 

iv ( 22 +1 \P) = \ JvU,RetrW^ 2 (x) 

eMS[U](w[f,w^)}. (22) 

As starting points for the iteration we chose the naive 
tree level values 

w$(0) = 1, u$G8) = 1. (23) 



In order to monitor the convergence of the above itera- 
tion, we introduce the error 

e («+D = \w& +1) (0) - I . (24) 

It turns out that measuring w^2 (0) is sufficient for mon- 
itoring convergence. In practice, the integrals in (|21I22|) 
are not calculated exactly. Only Monte-Carlo estimates 
w^\(3) with statistical errors (/?) are available. At 
the beginning of the iteration, it does not make sense 
to obtain a high precision estimate for an anyhow un- 
converged value of w\™\/3). We therefore adopted the 
following procedure: at the start of the iteration, only 10 
iterations are used to obtain the estimates wj" (0) and 

their statistical errors cr}™ (0). As soon as the error of 
convergence reaches the order of the statistical error, i.e. 

c<"> « ^(0) , (25) 

the number of iterations which are used for the estima- 
tors is increased by 10. The iteration stops when e- n > 
(and therefore also 0^(0)) drops below a certain num- 
ber which specifies the precision to be achieved for the 
parameters. Figure Q] shows the "thermalisation history" 
of the parameters w$ (0) /N c and (0) /N c as a func- 
tion of the total number lattice sweeps performed. Data 
are shown for (3 = 1.25 (SU(2)) and = 3.10 (SU(3)), 
which will turn out to correspond to a rather coarse lat- 
tice, and for = 1.55 (SU(2)) and = 3.5 (SU(3)), 
which arc in the scaling regime. After an initial oscil- 
lation, the estimators rapidly converge. Note that the 
spacing between two data points in figure [1] shows the 
number of lattice sweeps which were needed to estimate 
the integrals (|21|22|) . 
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FIG. 2: Static quark anti-quark potential using the present improved action for SU(2) (left) and SU(3) (right). 



In particular for small values of /3, several solutions of the 
non-linear equations (|19I20[) might exist. If one chooses 
to perform simulations in this regime of parameter space 

B. Static quark potential 

To investigate scaling, we will express the lattice spacing 
in units of the string tension a in order to calculate <x(/3). 
The static quark potential V{r) is so obtained from pla- 
nar Wilson loops. These loops are of size r x t, and the 
spatial links have been smeared to enhance the overlap 
with the quark anti-quark ground state. For the smear- 
ing procedure, we consider the spatial hypercube for a 
given time t: spatial links belonging to this cube are then 
cooled with respect to the 3-dimensional Wilson action. 
Cooling is performed by visiting each link of the lattice 
and replacing it the (normalised) sum of the adjacent 
staples. Instead of the Wilson action, one could use the 
3-dimensional version of the present action rather than 
Wilson's action. It turns out, however, that this choice is 
more time-consuming and does not produce better over- 
laps. 10 cooling sweeps through the lattice are performed 
to obtain one set of smeared links. Time like links are 
unaffected by smearing. The advantages of this smearing 
procedure are that it is easy to implement, it is fast com- 
pared with other smearing techniques and it is known to 
yield excellent overlap for SU(2) and SU(3) gauge theo- 
ries [13 and even for more exotic gauge theories such as 
for G(2) 0. 

In practice, the Wilson loops are fitted to a straight line: 

- In ^Re trW£**\ = V(r) t + const. , (26) 

where only data with t > t\ ow are included. This sup- 
presses the contribution from excited states. Because of 



the overlap enhancement, choosing t\ ow = 2a is sufficient: 
the linear t-fit represents the data with a x 2 /dof ss 2 or 
better for the /3 ranges explored in the present paper. We 
checked that larger values £i ow yield the same potential. 

In order to explore rotational symmetry breaking effects 
by the underlying lattice, "off-axis" distances for the 
quark anti-quark pair are considered as well. Potentials 
corresponding to cristalographic directions 

(100) (on-axis) , (110) (111) 

are taken into account. Our final result for the static 
potential is shown in figure [5] A fit of the on-axis data 
to 

o 

V{r) = V + g r (27) 

r 

is shown as well. The results of the fit for SU(2) and 
SU(3) are summarised in the table below: 





(3 


^Vconf 


Vo a a a a 2 


X 2 /dof 


SU(2) 
SU(3) 


1.45 
3.30 


400 
800 


0.527(2) 0.267(1) 0.0695(6) 
0.631(1) 0.317(1) 0.0666(3) 


1.1 
1.2 



Here, N con { denotes the number of independent lattice 
configurations used to estimate the Wilson loop expecta- 
tion values. Priority has been given to the SU(3) simula- 
tions because of their relevance for QCD. Because we are 
using an improved action with very good rotational sym- 
metry (see section W}. the point r = a can be included in 
the potential fit without hampering the value for x 2 /dof. 
A lattice sweep consists out of a Cabbibo-Marinari up- 
date followed by 4 reflections (for SU(2)) and 5 (for 
SU(3)), respectively. Each reflection replaces the actual 
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SU(2) 


P 




SU(3) 


P 


era 2 




2.20 


0.28(1) 




5.60 


0.278(6) 




2.25 


0.194(4) 




5.65 


0.219(4) 




2.30 


0.145(2) 




5.70 


0.171(2) 




2.35 


0.1022(9) 




5.75 


0.133(1) 




2.40 


0.0738(5) 




5.80 


0.1051(7) 




2.45 


0.0523(3) 




5.85 


0.0854(5) 




2.50 


0.0390(2) 




5.90 


0.0731(4) 




2.55 


0.0281(2) 




5.95 


0.0601(3) 




2.60 


0.0211(2) 




6.00 


0.0517(2) 










6.05 


0.0447(2) 










6.10 


0.0387(2) 



SU(2) Wilson action 



TABLE I: Measured scaling function a a (/?) for a 16 4 lattice 
using the Wilson action for SU(2) and SU(3) gauge theories. 



configuration by another one which possess the same ac- 
tion. We observe that this process enhances the ergod- 
icity of the algorithm: auto-correlations are reduced and 
a speed-up of thcrmalisation is observed. 20 "dummy" 
lattice sweeps are performed until the configuration is 
used for measurements. Especially for small (3 values, 
a smaller number of dummy sweeps might be sufficient. 
There is room for a further fine-tuning of the algorithm. 
Note that these sweeps are carried out for fixed Ki 2 i n 
(|12[) the values of which were determined during the ini- 
tial stage of thermalisation. 



IV. ASYMPTOTIC SCALING WITH 
IMPROVED ACTIONS 

A. The Wilson action - a case study 

The Wilson action has been widely studied and is widely 
used even nowadays. It has been known, however, for a 
long time that large deviations of the lattice spacing a(/3) 
from the perturbative scaling (|TJ) are common with this 
type of action. The purpose of the present subsection is 
to quantify this statement. 

The partition function employing Wilson action is given 
by 



VU^ exp{S wil [U]} (28) 
S wil [U] = (3 J2 ^-RetrW^ 1 ^). (29) 



Using the techniques outlined in the previous section, 
we calculated the static quark potential and the scaling 
function aa 2 ((3). The results, obtained from iV C onf = 800 
independent configurations on a 16 4 lattice, are sum- 
marised in table fl] Finite size effects are expected at 
the 1% level if the side length of the lattice exceeds 
1.5 fm 19]. For a 16 4 lattice, finite size effects therefore 
play a minor role as long as a a 2 > 0.04. 
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0.015625 
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2-loop 
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SU(3) Wilson action 
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FIG. 3: Scaling function aa (13) for a 16 4 lattice using Wilson 
action for SU(2) (top) and for SU(3) (bottom) gauge theories. 
Asymptotic scaling according to (|30I31|I is shown as well. 



Figure [3J visualises the data of table [T] In order to bring 
out any onset of asymptotic scaling, these data are com- 
pared with the perturbative scaling function at 1-loop 
and 2-loop level (see ([TT2|» : 



In [aa 2 ((3)] 
In [aa 2 (/?)] 



1-loop _ 4tt 2 _ 
asym ~ \P 



(30) 



2— loop 
asym 



In [aa 2 (&*)] 
Po 



(31) 



Mil JL 

P 2 Aef 



In 



(ftcf)] 



The perturbative scaling functions are normalised to re- 
produce the measured data for [3 — /3 re f- In figure [31 
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the highest considered value j3 is chosen for f3 rc {. It is 
remarkable that both for SU(2) and SU(3), the 2-loop 
scaling function ()3ip does not yield an improvement on 
the agreement of the data with the 1-loop formula (|3"0|) . 



B. Improved action 

Employing the procedure discussed in subsection IIII A| 
we generated well "thermalized" configurations (as well 
as the simulation parameters u>n(/3) and w 22 (/?)), see 
(fT5|) . for a range of (3 values which give reasonably sized 
lattice spacings for the present 16 4 lattice. The simu- 
lation parameters as well as the calculated value of the 
lattice spacing a in units of the string tension a are sum- 
marised in table 2 for the SU(2) gauge theory and in 
table 3 for the SU(3) case. The calculated scaling func- 
tions aa 2 (P) are shown in figure 4 and figure 5, respec- 
tively. As with the Wilson action, a comparison with 
the asymptotic scaling functions (|30|31| is made. An 
satisfactory agreement with asymptotic scaling on coarse 
lattices with aa 2 as large as era 2 w 0.1 is observed for 
both gauge groups. In both cases, the agreement with 
the 2-loop formula seems to be better than with the 1- 
loop result for aa 2 < 0.06. 



SU(2) improved action 




I 1 I 1 I 1 I 1 L_ 

1.2 1.3 1.4 1.5 1.6 
P 



Figure 4: Scaling function aa 2 (f3) for a 16 4 lattice us- 
ing the present improved action for SU(2). Asymptotic 
scaling according to (|30|31|) is shown as well. 







wu(P)/2 w 22 ((3)/2 N cont aa 2 



1.250 
1.275 
1.300 
1.325 
1.350 
1.375 
1.400 
1.425 
1.450 
1.475 
1.500 
1.525 
1.550 
1.575 
1.600 



0.62455(3 
0.63416(3 
0.64197(3 
0.64885(3 
0.65526(3 
0.66123(3 
0.66690(3 
0.67231(3 
0.67744(3 
0.68229(3 
0.68697(2 
0.69141(2 
0.69572(2 
0.69988(2 
0.70388(2 



0.13499(6 
0.15944(6 
0.17723(6 
0.19211(6 
0.20584(6 
0.21818(6 
0.22986(6 
0.24121(6 
0.25169(6 
0.26143(6 
0.27076(6 
0.27964(6 
0.28809(6 
0.29638(6 
0.30427(6 



400 
400 
400 
400 
400 
400 
400 
400 
400 
400 
400 
800 
800 
800 
800 



0.279(2) 
0.215(5) 
0.175(3) 
0.1473(7 
0.1244(7 
0.1068(9 
0.0922(7 
0.0787(5 
0.0695(6 
0.0599(3 
0.0528(3 
0.0452(8 
0.0400(3 
0.0351(4 
0.0311(2 







Wn(0)/3 w 22 (l3)/3 JV C o„f oa 2 



2.90 
3.00 



.10 
.15 
.20 

.25 
.30 
.35 
3.40 
3.45 
3.50 



0.58567(2 
0.60135(2 
0.61378(2 
0.61923(2 
0.62450(2 
0.62950(2 
0.63429(2 
0.63883(2 
0.64330(2 
0.64764(2 
0.65173(2 



0.11424(3 
0.14553(3 
0.16726(3 
0.17667(3 
0.18560(3 
0.19403(3 
0.20215(3 
0.20959(3 
0.21711(3 
0.22432(3 
0.23108(3 



800 
800 
800 
800 
800 
800 
800 
800 
800 
800 
800 



0.231(4) 

0.151(2) 

0.1122(7) 

0.0985(5) 

0.0851(4) 

0.0765(4) 

0.0666(3) 

0.0589(3) 

0.0532(2) 

0.0482(2) 

0.0424(2) 



Table 2: Simulation parameters of the improved action 
(fT7|) and the calculated scaling function <ra 2 (/3); SU(2) 
gauge theory, 16 4 lattice, AT con f independent configura- 
tions. 



Table 3: Simulation parameters of the improved action 
(fT7|) and the calculated scaling function aa 2 ((3); SU(3) 
gauge theory, 16 4 lattice, N con f independent configura- 
tions. 



8 



SU(3) improved action 



0.25 



0.125 



0.0625 



0.03125 




_L 



2.8 2.9 



3.1 



3.2 
P 



3.3 3.4 3.5 3.6 



Figure 5: Scaling function aa 2 {(i) for a 16 4 lattice us- 
ing the present improved action for SU(3). Asymptotic 
scaling according to (|3QI3ip is shown as well. 



C. Comparison with other actions 

In this subsection, the more important case of a SU(3) 
gauge group is investigated. Two popular actions which 
do not involve tadpole improvement, but invoke a renor- 
malisation group investigation, are the RG-Iwasaki ac- 
tion @, @| and the DBW2 These actions are of the 
type 

S[U] = P Yl [co^RetvW^ix) 



Cl ± RetrW ^ {x) 



(32) 



and differ by the choice of ci (Note that cq = 1 — 8ci for 
a proper definition of the bare gauge coupling) : 



ci 
ci 



-0.331 
-1.4088 



(RG-Iwasaki) 
(DBW2). 



A detailed investigation of the scaling properties of these 
actions can be found in [2(J. We will here focus on their 
properties concerning asymptotic scaling. 
We will need the lattice spacing a in units of the string 
tension. For the case of the RG-Iwasaki action and the 
DBW2 action, data for a/ro with the Sommcr parameter 
ro are taken from the work by Necco [2(J. Using r$ w 
0.5 fm and yfo w 440 MeV, a factor 

ar\ k, 1.21 

is used to convert a 2 /r§ to a 2 a. 

In order to study whether the present improved action 
(see equations (|16H18p ) is superior to an action with stan- 
dard tadpole improvement, we here also study the "2 x 2" 



action with tree-level coefficients and standard tadpole 
removal: 



S[U]( W11 ) = p [ 



l 



48U&G8) 



3 w n 09) 
RetvW% 2 (x) 



RetrW** 1 ^) 



(33) 



where w\\ must be self-consistently calculated from 

w n (l3) = ^ VU^RetrW^ix) exp{S[U](w u )} . 

_ (34) 
This action was used in [2 11 ] to study thermodynamics. 
There it was observed that tadpole improvement largely 
reduces the cutoff effects which hamper the calculation of 
the pressure and the thermal energy density in the SU(3) 
high temperature phase. Our findings for w\\{(3) and for 
the scaling function aa 2 (/?) are summarised in table IIVI 



p 


3 ton 09) 


iVconf 


aa 2 


2.60 


0.57313(3) 


600 


0.115(1) 


2.70 


0.58559(2) 


600 


0.0875(6) 


2.80 


0.59715(2) 


600 


0.0683(4) 


2.90 


0.60783(2) 


600 


0.0520(3) 


3.00 


0.61784(2) 


600 


0.0431(3) 



TABLE IV: Simulation parameter of the standard tadpole im- 
proved action (133[) and the calculated scaling function aa?{(5); 
SU(3) gauge theory, 16 4 lattice, N con { independent configu- 
rations. 



For a more quantitative investigation of asymptotic scal- 
ing, the deviation from asymptotic scaling is measured 
by the ratio: 



R(P) 



2 W 



(35) 



where the lattice spacing squared a 2 is either provided in 
units of the string tension or in units of the Sommer pa- 



). The function 



asyni 



311. Because the defini- 



rameter ro (as e.g. done in [2C 
is provided at 2-loop level by 
tion of a 2 sym (P) involves an arbitrary normalisation, the 
absolute value of R in (|35f is meaningless. The data are 
normalised such that R = 1 is attained for the smallest 
lattice spacing considered. Asymptotic scaling will from 
([33)1 be signalled by the function R(f3) becoming flat for 
sufficiently large values of j3. Since the absolute size of 
the bare gauge coupling g (and therefore of = Q/g 2 ) 
depends on the details of the regularisation scheme and 
the action, we will study R as a function of the lattice 
spacing squared in physical units. 

The results are shown in hgure 3] for the actions un- 
der investigations. It turns out that any sort of tadpole 
improvement largely improves scaling along the lines of 
the asymptotic formula. The Wilson action, but also 
the DBW2 and Iwasaki actions, show large deviations 
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FIG. 4: Approaching asymptotic scaling using several ac- 
tions. 



from asymptotic scaling. In contrast, both, the standard 
and our new tadpole improved action, seem to perform 
equally well. A possible explanation could be that both 
actions remove the Oq irrelevant terms such that asymp- 
totic scaling sets in for rather coarse lattice spacings. The 
next section will, however, reveal that this interpretation 
is only justified for the new action proposed in this paper. 



V. ROTATIONAL SYMMETRY BREAKING 

The irrelevant terms of Oq in ((5J) are built up from ex- 
pressions which explicitly violate rotational symmetry. 
Hence, the absence of the Oq terms can be checked by 
calculating the amount (for a proper definition see be- 
low) of rotational symmetry breaking as a function of 
the lattice spacing a. 

The present improvement scheme belongs to the class 
of tree level improvements: it relies on the expansion 
(JTUJl of the action in powers of the lattice spacing a. A 
cancellation of the Oq terms can be hampered by terms 
which depend logarithmically on a. The cancellation can 
still be made complete if loop corrections are considered 
as well as tadpole improvement. At the present stage, 
there are two crucial questions: Are loop corrections still 
large for the present range of lattice spacings? Is the new 
improvement scheme superior to the standard approach 
so that the additional complexity of the new scheme is 
justified? 

For an answer to these question, we need to quantify 
the amount of rotational symmetry breaking. For this 
purpose, we invoke the method introduced by the QCD 



FIG. 6: The measure 5 V of rotational symmetry breaking as 
a function of the scaling function aa 2 for SU(3) gauge theory 
with different actions. 



Taro collaboration in [la] . Let V on (r) denote the "on- 
axis" static quark potential obtained from quarks posi- 
tioned along the main crystallographic direction previ- 
ously called the (100) direction. Data for which r is not 
an integer multiple of a are made available by means of 
the fit (|27|) . Let furthermore call V(r) all data arising 
from quarks positioned along the (110) and (111) direc- 
tions. These data are called the "off- axis" data. 5V(r) 
denotes their statistical errors. With these definitions, 
the measure of rotational symmetry breaking is given by: 



£ 

off 



[V(r) - Von(r)} 2 
V(r) 2 SV 2 (r) 



/(£ 



i 



5V 2 {r 



■) , (36) 



where the sum extends over all "off-axis" data. 
Figure [5] shows S v as a function of the lattice spacing 
squared. For the case of the Wilson action, we observe a 
linear scaling of 5 V with a 2 : 



S v 



(Wilson action) 



In the case of the action employing standard tadpole 
improvement, the symmetry breaking effects are signifi- 
cantly reduced, but, still, S v rises linearly with a 2 . Using 
the new action proposed in this paper, we find an addi- 
tional drastic reduction of rotational symmetry breaking 
effects. Moreover, it seems that this time the functional 
dependence of 8 V on a 2 seems to change qualitatively: 
The data now indicate that S v is of higher order in a 2 : 



2 4 

oc a a 



(SU(3) imroved action) 



However, further numerical investigations are necessary 
to support this claim. If it is supported by numerical 
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simulations, it would imply that only the present action 
can completely remove the irrelevant 6 contributions. 
It is already clear that the standard tadpole improved 
action certainly fails this task. 

Note also that the data suggest that the Wilson and stan- 
dard tadpole data approach a small but finite value at 
era 2 = 0. The size of this value clearly depends on the 
amount of tadpole contributions. Our preliminary inter- 
pretation of this finding is as follows: contributions of 
tadpole loops solely arise in lattice regularisation and, 
therefore, add substantially to the amount of rotational 
symmetry breaking present in the static potential. In ad- 
dition, tadpole loops are generically UV divergent. This 
might lead to a small, but finite value of S v even for very 
small values of the lattice spacing. Further investigations 
are clearly needed to settle this question. 

VI. CONCLUSIONS 

The properties of improved actions with respect to 
asymptotic scaling has been thoroughly investigated in 
this paper. A focal point of the present study is tad- 
pole improved tree level actions. A new scheme for tad- 
pole improvement has been proposed and it has been 
contrasted to the heuristic tadpole approach, which is 
standard in the literature. It has been shown that the 
standard tadpole scheme is a mean field approximation 
to the scheme proposed here. 

The numerical results for the scaling function era 2 {(3) re- 



veal that both types of tadpole improved actions yield 
results of equal quality as far as asymptotic scaling is 
concerned. By contrast, loop improved actions (which 
do not make use of tadpole improvement) produce much 
bigger deviations from asymptotic scaling. 

The amount S v of rotational symmetry breaking (see 
(|36[) )in the static quark potential was used to compare 
the quality of both tadpole improvement schemes. Al- 
though the function S v (a) is much smaller for the stan- 
dard tadpole action than for the Wilson action, the func- 
tional dependence on the lattice spacing a is the same in 
both cases. In contrast, we have seen first numerical ev- 
idence that S v is of higher order in a 2 if the new tadpole 
improved action is used. We here argue that the generic 
tadpole scheme fails to eliminate the leading order irrel- 
evant terms of the action. The data indicate that the 
new action cancels these terms from the action for the 
range of lattice spacings considered without taking into 
account loop corrections. 

Our approach to tadpole improvement can in principle be 
extended to kill off the next to leading order irrelevant 
terms as well. The question whether loop corrections 
must be considered then is left to future work. 
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